The accuracy of roughness exponent measurement methods 
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We test methods for measuring and characterizing rough profiles with emphasis on measurements 
of the self-afBne roughness exponent, and describes a simple test to separate between roughness 
exponents originating from long range correlations in the sign signs of the profile, and roughness 
exponents originating from Levy distributions of jumps. Based on tests on profiles with known 
roughness exponents we find that the power spectrum density analysis and the averaged wavelet 
coeflicients method give the best estimates for roughness exponents in the range 0.1 to 0.9. The error- 
bars are found to be less than 0.03 for profile lengths larger than 256, and there are no systematic bias 
in the estimates. We present quantitative estimates of the error-bars and the systematic error and 
their dependence on the value of the roughness exponent and the profile length. We also quantify how 
power-law noise can modify the measured roughness exponent for measurement methods different 
from the power spectrum density analysis and the second order correlation function method. 

PACS numbers: 68.35.Ct,05.40.Fb 
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I. INTRODUCTION 



In 1990 Bouchaud et al. [l| proposed that the rough- 
ness exponent for profiles from three dimensional fracture 
surfaces was universal - independent of material proper- 
ties and fracture mode. This has since been the working 
hypothesis in the community studying fracture rough- 
ness. To claim universality one needs good measurements 
and to be aware of the inherent biases and limitations of 
the different methods used for measuring the roughness 
exponent. The use of several independent methods are 
a prerequisite for a good estimate of the roughness ex- 
ponents. This is even more true now as the study of 
self-afHne surfaces goes beyond simple measurements of 
the roughness exponent and the higher order statistics 
and the shape of the height-difference distribution func- 
tion p(A/i, I) are also studied. 0, 0] Now a thorough 
measurement of the self-affinity of a surface should also 
include checks for corrections to scaling and a survey of 
the higher order statistics i.e. multi-scaling, and if needed 
a check for anomalous scaling. 

Research on surfaces morphology with focus on the 
scaling properties of the surfaces has beenpursued since 
the work of Mandelbrot et al. in 1984. % These stud- 
ies have not been restricted to fracture surfaces. Exam- 
ples of surfaces that have been studied during the last 
twenty years are fluid fronts in disordered media, fire 
fronts in paper, atomic deposition surfaces, fracture sur- 
faces, DNA base-pair sequences and the time signal of 
the heart rhythm. These surfaces have been shown to 
have statistically self-affine scaling properties. The self- 
affine scaling is an anisotropic scaling of the system. This 
is seen in the scaling of the height-difference distribution 
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FIG. 1: (Color on-line) If a surface of a given linear size L\ is 
with 1/2 = AZ/i in the horizontal direction, the surface must 
be scaled with h{Xx) — X~''h{x) in the vertical direction to 
be statistically similar. is the roughness, or Hurst exponent, 
measuring the degree of anisotropy. C = 1 gives self-similar 
scaling. In this figure the part of the black profile which is 
inside the square is rescaled with the scaling relation above. 
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An example of self-affine scaling in two dimensions is 
shown in Fig. [TJ Systems with a common roughness ex- 
ponent are said to belong to the same universality class, 
and are therefore controlled by the same fundamental 
physical law. An early summary of different universality 
classes and the different surfaces studied can be found in 
Barabasi and Stanley . 

The surfaces described above are constructed in differ- 
ent ways. Restricting the discussion to two-dimensional 
surfaces or profiles we can divide the profiles into three 
groups. The first group of profiles grow from an initial 
planar profile or line into a rough profile. Examples from 
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this group of profiles are wetting fronts, deposition fronts 
and three dimensional fracture surfaces constrained to 
grow in a plane. The second group of profiles grow at 
the end points. Examples from this group of profiles are 
DNA base-pair sequences, any times series and fractures 
growing in two-dimensional systems. The third group 
is made by coalescence of micro cracks. Examples from 
this group of cracks is crack growth by plastic micro void 
coalescence. The initially flat profiles are often character- 
ized by the Family- Vicsek scaling relation of the profile 
width w{L,t) cx L"f{t/L^) 0], where a is the roughness 
exponent measured when the system is saturated and 
z = is the dynamic exponent. f3 is the growth expo- 
nent characterizing the growth of the profile width before 
saturation 'w{L,t) cx The values of the three exponents 
then give the universality class. Profiles that do not obey 
Family- Vicsek scaling are said to have anomalous scaling. 

The concept of anomalous scaling can not be applied 
to the second group of profiles as they are only 1 -f 1 
dimensional, while the profiles in the first group are 2 + 
1 dimensional. The fracture profile in two dimensions 
are strictly speaking not 1 + 1 dimensional since there 
are two space dimensions and no time dimension. Each 
single part of the fracture is frozen once it is placed in the 
system. It is therefore no time evolution of the fracture 
in 2 -f 1 dimensions. 

The aim of this paper is to give a survey of the dif- 
ferent methods used for measuring the roughness expo- 
nent and give an estimate of the expected errors and 
biases. This paper will use some of the methods used 
by Schmittbuhl et al. in excluding the fractal mea- 
surements and the return methods since these methods 
have already been tested there. We repeat the test done 
with the local window methods and the power spectrum 
density analysis. Other methods such as the detrended 
fluctuation analysis 8], the averaged wavelet coefhcient 
method and the height difference distribution are also 
included. To accurately asses the results from the dif- 
ferent methods the methods must be applied to profiles 
with known roughness exponent. Both the Voss method 
[Tot ^^'^ wavelet method [III for generating profiles 
have been used and analyzed. We will also study the 
systematic error introduced by adding a power law noise 
to the signal, and discuss that both long range correla- 
tions in the sign change and a power law noise can give 
surfaces with roughness exponents. This is motivated 
by recent results for the central force and fuse model of 
fracture. [13, [3 ■ Our interest in the phenomena of self- 
afhne surfaces comes from the studies of fracture surfaces 
both experimental and numerical, and the difhcultics we 
have encountered while measuring the roughness expo- 
nent. 

The structure of this paper is as follows: In Sec. Ullthe 
different measurement methods are presented. In Sec. lIIII 
we presents the results from the measurements done on 
the generated profiles. Samples of the power law scal- 
ing for each method is presented and compared with the 
other methods. The results for both the generation meth- 



ods described in the last paragraph are shown. We find 
that the power spectrum density analysis and the aver- 
aged wavelet coefficient method give the most accurate 
for roughness exponents in the range 0.1 — 0.9. In Sec. 
II VI we study the size dependence on the systematic errors 
for the different methods. We observe as expected that 
the errors decreases for larger system sizes, but several of 
the measuring methods have relatively large systematic 
error and the size of this error depends on the measured 
roughness exponent. In Sec. |V]we show that power law 
noise might disguise as self-affinity when profiles are ana- 
lyzed with some methods. We then describe how one can 
separate the contribution from the power law noise and 
the sign change correlation to the roughness exponent. 
And finally in Sec. I VII we give a summary and come with 
some recommendations to researchers studying surface 
roughness. 



II. MEASUREMENTS OF THE ROUGHNESS 
EXPONENT 

Until recently one used different methods for measur- 
ing the roughness exponent on experimental and numer- 
ically made fracture surfaces. On surfaces from experi- 
ments one applied local (or intrinsic) methods which use 
scaling with a (intrinsic) length scale I much smaller than 
the system size L or spectral methods like the power 
spectrum density. On the other hand, for surfaces from 
numerical simulations one applied global (or extrinsic) 
methods which use scaling with the system size. This 
was was done by necessity as in experiments it is not easy 
to make samples of the same material with sizes ranging 
over several orders of magnitude, or experimental setups 
that can do measurements over the same orders of mag- 
nitude. While for numerical simulations the restriction in 
computing power made it difficult to create large enough 
samples to use the local methods. But during the last few 
years and due to the growth in computing power numer- 
ical simulations of fracture have produced large enough 
samples in such numbers that both the real space and 
spectral local methods are used in numerical studies of 
fracture surfaces. 

Below we will describe some local and global methods 
which are used today for measuring the roughness expo- 
nent. All the methods we consider in this paper work 
on profiles. I. e. profiles in two-dimensional planes cut 
from a surface embedded in three dimensions or a one- 
dimensional trace or path embedded in two-dimensions. 

The local window methods all measure the scaling of 
a characteristic width as a function of the window size. 
The width is a measure of how large the fluctuations of 
the profiles in a window of length I are. We will look 
at three different methods which use different definitions 
of the characteristic width: 1. The variable bandwidth 
method (VB) 

wvB{l)^i{ihix)-h)^)L)'/^ o.l'^, (2) 
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where h = {h{x)). 2. The detrended fluctuation analysis 
(DFA), where the local linear trend in each window is 
subtracted 



WDFA{l)^i{{h'{x)^h'f)L)'/^ 0,1'^. 

3. the maximum - minimum method (MM) 



(3) 



WMAiil) = {max{h{x),x e {xo,Xo + 1} 

—min{h{x),x&{xo,Xo + l})LCcl''. (4) 

means averaging over the system size L. Aver- 
aging over different samples is implied. The roughness 
exponent C is then found as the power law scaling with I 
of the characteristic width. 

An alternative measuring method is to look at the sec- 
ond order correlation function 

C2{r)^{{h{x~r)-h{xyf)]!' ^r'^; (5) 

The roughness exponent can also be found from the 
two spectral methods as the scaling of the power spec- 
trum P{k), [ij] which is the Fourier transform of the 
correlation function Eq. with the wave number for 
the power spectrum analysis 



P(fc)afc-2(';+i), 



(6) 



and as the scaling of the averaged wavelets coefficients 
MK[/i](a) as a function of the scale variable a 



W[h]{a) cx a 



+C+1/2 



(7) 



The global average width method is similar to the in- 
trinsic one, but with the window size I substituted with 
the system size L. 



W{L) ^ {h{x) - /i) cx L^, 



(8) 



Similarly we have the global maximum-minimum method 

Wmm{L) = {max{h{x),x e {0,L} 

-mm{h{x),x e {0, L}) cx L^. (9) 

The roughness exponent measured by the global meth- 
ods give the same results as the local methods only of the 
scaling of W{L) is the same as the scaling of w{l), I ^ L. 
One example of when this is not the case is when the 
profiles shows anomalous scaling, In this paper we 
will not consider the global methods. 

In two recent papers the question of whether the frac- 
ture surfaces measured in experiments are self-afhne or 
multi-afhne has been discussed. 0,11] In Santucci et al. Q 
two slightly different methods for measuring the fracture 
roughness have been discussed. The basis for these meth- 
ods is to assume that the distribution function p{Ah, I), 
Ah{l) = h{x + I) — h{x) is 'Gaussian-like'. The first 
method was introduced in the studies of directed poly- 
mers in random media, p^ . The fc-th moment of h{x). 
is defined below 



Ck{i)^{\h{x + i)-h{xty^'' 



To check whether or not the surface have self-affine or 
multi-affine scaling one calculates the ratio of the fc-th to 
the second moment 



Rk{l) 



{\h{x + i)~hix)\''y/'' 



(11) 



({h{x + 1) - h{x)yy/^ 

which for a true Gaussian distribution should reduce to 



i/fc 



(12) 



By plotting Rk/R'^ one should then get straight lines 
if the surface is self-affine since Cfc = when the sur- 
face is self-affine. If these lines should fall on top of each 
other, the underlying distribution is also Gaussian. The 
roughness exponent can be found by finding the slope of 
Cfc/C^, where is the fc-th moment for a true Gaus- 
sian distribution, in a double logarithmic plot, which will 
return the same roughness exponent as C2{r). 

The other method is to construct the distributions 
p{Ah, I) of the height difference over distance I and plot- 
ting p{Ah,l)V2na'^ versus {Ah)/\/2a^^ where is the 
fluctuations of Ah over I. For Gaussian distributions the 
plot should be a parabola pointing downward in a semi- 
logarithmic plot. The roughness exponent can then be 
found from the power law scaling of slope of the fluctua- 
tions in Ah, a{l) cx l'^. 

Many of the methods described above involve the scal- 
ing of a characteristic width with a window length. At 
least two mechanisms can lead to the increase in the char- 
acteristic width. [31 The flrst one is that there are spatial 
correlations in the sign change of the steps in the profile, 
and the second one is a Levy like jump distribution. To 
check how these two mechanisms contribute to the mea- 
sured effective roughness exponent on can do two differ- 
ent modifications to the measured profiles. If one sets 
the size of each jump equal to unity using 



ho{x) 



lim 



sgn{hix))\h{x)\idx, 



(13) 



one will only measure the characteristic width that are 
caused by the correlations in the sign changes as any in- 
formation carried by the amplitude will be removed. It 
has been shown [l7] that for ho{x) the roughness expo- 
nent is 



C = max(i,C/i) 



(14) 



(10) 



where (^h is the roughness exponent of h{x). Thus for 
a roughness exponent less 1/2 some information in the 
roughness exponent will be in the jump distribution. 

To remove any correlation that might be in the sign 
changes, but keep the information in the jump distribu- 
tion one can randomly rearrange the position for each 
jump. The characteristic width measured on this ran- 
domly rearrange profiles hr{x) will now only depend on 
the jump distribution, and therefore the measured rough- 
ness exponent is due to the jumps. 
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III. ACCURACY OF THE MEASUREMENT 
METHODS 

In Schmittbuhl et al. Q the reliability of self-affinc 
measurements is addressed for some of the methods 
above. We will repeat the measurements for the power 
spectrum density analysis and local window methods, 
and compare with the additional methods from Sec. |TT] 
to check the validity of the roughness exponents we mea- 
sure from the generated self-afHne fracture surfaces. 

We base our choice of good measurement methods on 
two criteria: 1. The fitting of the power law should be 
good i.e. good linearity in the double logarithmic plots. 
2. The systematic error for the method should be small 
and stable over a range of different roughnesses. 

We generate artificial surfaces u sing two methods. The 
first method is the Voss method and the second is 
the wavelet method 'Tl'|. The profiles are generated by 
first creating long profiles of a length much larger than 
the profile sizes we want to study. We then cut out a 
piece of desired length at a random position from the 
long profiles. In this study the long profiles had a length 
of 65536. Before we measure the roughness exponent we 
remove any linear trend in profiles. The tests are done 
for C G {0.1,0.9} and system sizes L < 16384. For each 
roughness exponent and each system size 100 samples 
were generated. Samples of these profiles are shown in 
Fig. [21 Before one starts to measure the roughness ex- 
ponent it is from our experience always a good idea to 
visually inspect the profiles and compare them to self- 
affine profiles with known roughness exponents. 



In Figs. [3] to [7] we present the different measurements for 
L = 512. This system size is now numerically accessi- 
ble for many different numerical models of fracture and 
large enough to give good estimates for the roughness ex- 
ponent using the local methods. The profiles were made 
with C = 0.6 and all straight lines in the figures repre- 
sents this C" value. These plots arc from the Voss profiles. 
From Fig. [3] one see that for the local window methods 
in Eq. ^ to Eq. (jH) the detrended fluctuation analysis 
gives the best estimate of the roughness exponent. The 
detrended fluctuation analysis gives a straight line in the 
double logarithmic plot from I = 16 to I = 512 while 
the variable bandwidth and the Max-Min data are more 
curved compared to the C — 0.6 lines. 
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FIG. 2: (Color on-line) Sample profiles from the Voss algo- 
rithm (a) and the wavelet algorithm (b) . The profiles are from 
top to bottom for roughnesses 0.2, 0.5 and 0.8. To make the 
scale of the fluctuations the same for samples with different Cj" 
the samples are rescaled with W(;{L)/Wo.5{L). 
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The power spectrum density method in Fig. |4] as well 
as the averaged wavelet coefficients method in Fig. [5] 
both give good estimates of the roughness exponent. For 
Ck/C^ we can see in Fig. IHlthat the second order correla- 
tion function gives an estimate of the roughness exponent 
comparable which is below the correct value. In addition 
we see that the collapse of the different moments is show- 
ing that the profiles are self-affine, not multi-affine. Fig. 



FIG. 3: (Color on-line) Measurements of the roughness expo- 
nent using intrinsic window methods. From the top: The local 
MAX-MIN method Eq. (g]), the variable bandwidth method 
Eq. ([2]) and the detrended fluctuation analysis Eq. [S] 



[7] shows the scaling of a{l) which give results similar to 
the second order correlation function. 



5 




100 p 




1000 



FIG. 4: (Color on-line) Measurements of the roughness expo- 
nent using the power spectrum density analysis, Eq. (6) 
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FIG. 5: (Color on-line) Measurements of the roughness expo- 
nent using the averaged wavelet coefficients method Eq. [T] 
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FIG. 6: (Color on-line) Measurements of the roughness expo- 
nent using second order correlation function Eq. ([5]). 



FIG. 7: (Color on-line) Measurements of the roughness expo- 
nent using the scaling of a{l) for the distribution p{AH,l) 



the wavelet method in Fig. [S] we notice some differ- 
ences. Using the power spectrum density analysis and the 
averaged wavelet coefficients methods the profiles made 
by both methods we obtain good measurements of the 
roughness exponent. The results for when using the av- 
eraged wavelet coefficient method on the wavelet gener- 
ated profiles were no surprise as this should a priori give 
correct results. When one use the detrended fluctuation 
analysis and the variable bandwidth, the Max-Min and 
C2 methods we see a more pronounced deviation from 
the true roughness exponent used to create the profiles 
for the highest C for the wavelet generated profiles. 



This result suggests that the profiles made by the wavelet 
method is not as good as the Voss methods for generat- 
ing self-affine profiles. A test of the self-affinity of the 
profiles confirms this. In Fig. [10] we clearly see that the 
wavelet generated profiles have corrections to scaling at 
small scales. This correction to scaling can also be seen 
for p{Ah, I) in Fig. [TT] For the high ^-values the wavelet 
generated profiles the tail is broader than a Gaussian 
distribution. This effect becomes larger for the largest 
Cs. The profiles generated by the Voss algorithm is more 
narrow than the Gaussian distribution, and also for these 
profiles is the effect greatest for the largest values. In 
the rest of this paper we will only consider the profiles 
made by the Voss method. 



When we compared the result for the profiles made with 
the Voss algorithm in Fig. [5] and the profiles made with 



From Figs. [8] and [9] one can see that the power spec- 
trum density analysis and the averaged wavelet coeffi- 
cients method give the best estimates for the roughness 
exponent in the range 0.1 — 0.9. For these two methods 
the test measurements done here gave error bars from 
the power law fits equal or less than 0.03. The detrended 
fiuctuation analysis gave good estimates above C = 0.5, 
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FIG. 8: (Color on-line) (measured — (true for profiles generated 
by the Voss algorithm, a) The detrended fluctuation analysis 
(♦), the variable bandwidth (•), the Max-Min (■) and the 
standard deviation of Ah{l) (+). b) Second order correla- 
tion function (♦), power spectrum density analysis (•) and 
averaged wavelet coefficients (■). 
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FIG. 9: (Color on-line) (measured — (true for profiles gener- 
ated by the wavelet algorithm, a) The detrended fluctuation 
analysis (♦), the variable bandwidth (•), the Max-Min (■) 
and the standard deviation of Ah{l) {+). b) Second order 
correlation function (♦), power spectrum density analysis (•) 
and averaged wavelet coefficients (■). 



and systematically to high estimates for C < 0.5. The 
other local methods did only give good estimates around 
C = 0.5, and showed an systematic drift towards lower 
C values for ( > 0.5, and towards higher ( values for 
( < 0.5. Of the global W(L) gave the most accurate 
estimates, but gave systematically to high estimates for 
C < 0.5. 

The results for the roughness exponent were done with 
the following limits on for the regions where we did a 
least square fit. For the variable bandwidth method, the 
Max-Min method and detrended fluctuation analysis the 
regression was done in the range / € [16, L/2]. For the 
second order correlation function the regression was done 
for I e [2,L/8]. For the power spectrum analysis the 
regression was done for k G [0,0.2]. For the averaged 
wavelet coefficients method the regression was done for 
a G [8, L/4]. For the standard deviation of the height- 
difference the regression was done for I G [1, L/16]. 

We conclude this section by checking that the rough- 
ness exponent we measure from the profiles generated 
by the Voss algorithm that the self-afFine scaling origi- 



nates from the long range correlations in the sign changes. 
From the original profiles h{x) we construct the profile 
ho{x) where all jumps are on equal size. When ( > 0.5 
one sees from Fig. [T^ the same general results as pre- 
sented above for h{x). The power spectrum analysis, the 
averaged wavelet coefficients method and the detrended 
fluctuation analysis give the best estimates. However all 
estimates are from these three methods are now system- 
atically 0.05 to low compared to the true value. Also the 
global methods give good estimates, but these are also to 
low. The other local methods show the same large errors 
as they did for h{x). 

When C < 0.5 the picture is is not that clear. The re- 
sults for all the different methods show that they measure 
values above (h, but below 1/2 showing that not all the 
self-affine information for < 0.5 is in the sign change 
correlations for the self-affine profiles generated with the 
Voss algorithm. 
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FIG. 10: (Color on-line) Rk/Rk for profiles of length 512 
made by a) the Voss algorithm and b) the wavelet method for 
C = 0.7. 
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FIG. 11: (Color on-line) p{Ah,l) for profiles of length 512 
made by a) the Voss algorithm and b) the wavelet method 
for ^ = 0.7 Three different length scales are presented here: 
I = 64 (♦), / = 128 (•) and / = 256 (■). 



IV. SIZE DEPENDENCE 

In the last section we described how the different mea- 
surement methods behaved for profiles with L = 512. In 
this section we will expand our study by considering how 
the roughness exponent estimates change with the sys- 
tem size. The accuracy with which one can measure the 
roughness exponent is dependent of the length of profiles. 
In [7| Schmittbuhl et al. reported on the size dependence 
of the error bars for the variable bandwidth method, the 
MAX-MIN method and the power spectrum analysis. We 
have performed the measurements of the roughness expo- 
nent with the methods described above on samples gen- 
erated by the Voss algorithm for different system sizes. 
The results are presented in Fig. [13] to Fig. [2TJ 

As the system size increase the error decreases for the 
all the measurement methods. But even for large system 
sizes the different methods measure the with systematic 
errors which are dependent on the value of the C- These 
systematic errors are similar to the ones shown in Fig. [S] 

The detrended fluctuation analysis, the power spec- 
trum density analysis and the averaged wavelet coeffi- 
cient method all have systematic errors smaller that 0.05 



for L > 256 for all values of C. The local window methods 
and the second order correlation function method over- 
estimate C for C < 0.5, and under-estimate C for > 0.5 
as reported in Sec. IIIII 
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FIG. 12: (Color on-line) ^measured — Cfue for ho{x) for pro- 
files generated with the Voss algorithm. L = 512. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of A/i(i) 
b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefficients 
(■)• 
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FIG. 13: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for C,true = 0.1. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of A/i(/) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefficients 
(■)• 



I 




00 CO (M 
^ 

o o ^ 

(M ^ 00 



00 
00 
CO 



(a) 



I 




CO 



00 


CO 


(M 




00 


CO 










i-H 








05 


00 


1-H 


(M 


LO 


o 


o 


o 


i-H 


CO 








1 — 1 






00 


CO 



(b) 



FIG. 14: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.2. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 
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FIG. 15: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.3. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 



FIG. 16: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.4. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 
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FIG. 17: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.5. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 



FIG. 18: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.6. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 
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FIG. 19: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.7. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 



FIG. 20: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.8. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 
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FIG. 21: (Color on-line) Systematic errors in roughness expo- 
nent measurements versus system size for (true = 0.9. a) The 
detrended fluctuation analysis (♦), the variable bandwidth 
(•), the Max-Min (■) and the standard deviation of Ah{l) 
(x). b) Second order correlation function (♦), power spec- 
trum density analysis (•) and averaged wavelet coefflcients 
(■)• 
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FIG. 22: (Color on-line) Two samples profiles of length 1024, 
the one with solid lines are made with the Voss algorithm 
with (^true = 0.7, and the dashed one is a random walk with 
Levy-distributed jumps with a — 1.5 



V. POWER LAW DISTRIBUTED STEP SIZES 

We will in this section compare two different sets of 
profiles which have the same roughness exponent when 
measured with the detrended fluctuation analysis, but 
different when measured with the power spectrum den- 
sity analysis. In Fig. [221 one sample profile from each 
set of profiles are shown. One profile made with the Voss 
algorithm with Ctrue = 0.7, and one profile which is a ran- 
dom walk with Levy-distributed jumps with a = 1.5 are 
shown. As seen in Fig. [22] the main difference between 
these two profiles are the existence of the large jumps 
in the second profile. When one use the detrended fluc- 
tuation analysis to measure the roughness exponent for 
these different sets of profiles it gives the same roughness 
exponent for both sets as shown in Fig. [53) The only 
difference is that the amplitude of wdfa{1) is larger for 
the Levy-flight profiles than for the Voss profiles. If one 
use the power spectrum analysis, one will measure two 
different roughness exponents. ( — 0.7 for the proflles 
made with the Voss algorithm, and ( = 0.5 for the Levy 
profiles. By using the two modified proflles ho{x) and 
hr{x) described in Sec. [Til and measuring the roughness 
exponent on these one can distinguish between the two 
different sets of profiles. From Table [J one sees that the 
power spectrum density analysis do not measure C 0.5 
for the random walk with Levy-distributed jumps. De- 
trended fluctuation analysis will measure — 0.7 except 
when the jumps are all set to the same size. For the 
profiles generated with the Voss algorithm both the de- 
trended fluctuation analysis and the power spectrum den- 
sity analysis give C = 0.7 except for when the long range 
correlations are destroyed in hr{x). This again show that 
the roughness exponent for C, > 0.5 for the Voss gen- 
erated proflles comes from the correlations in the sign 
change. It also shows that the effective roughness expo- 
nent measured with the detrended fluctuation analysis on 
the Levy proflles are caused by the power law noise and 
not the sign change correlation. In addition one notes 



FIG. 23: (Color on-line) Roughness exponent measured with 
detrended fluctuation analysis for a proflle made with the Voss 
algorithm with C^true ~ 0.7(top) and from a random walk 
with Levy-distributed jumps with a — 1.5(bottom).Both the 
straight dashed lines in the figure show = 0.7. 100 samples 
of length for each type of profiles were used in the calculations 
of the roughness exponent. 
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FIG. 24: (Color on-line) Roughness exponent measured with 
PSD for a proflle made with the Voss algorithm with (^tme = 
0.7(bottom) and from a random walk with Levy-distributed 
jumps with a = 1.5(top). The corresponding straight lines 
are for C, = 0.5 (top) and C, = 0.7 (bottom). 100 samples of 
length for each type of profiles were used in the calculations 
of the roughness exponent. 



that different measurement methods give completely dif- 
ferent roughness exponents in this case. 



TABLE I: Table of roughness exponent measured with the de- 
trended fluctuation analysis and the power spectrum analysis 
on the original proflles (("), the modifled profiles with equal 
jump size (Co) and the randomly rearranged proflles (Cr). 



Method 



DFA 



PSD 



Proflle set C Co Cr C Co Cr 
Voss 0.7 0.7 0.5 0.7 0.7 0.5 
Levy 0.7 0.5 0.7 0.5 0.5 0.5 



The power law noise give different corrections to the 
different measuring methods describe in Sec. [TT] as was 
shown earlier in this section. To quantify these correc- 



15 




0.5 1 1.5 2 2.5 3 3.5 



a 

FIG. 25: (Color on-line) Roughness measurements on random 
walks with power law distributed jumps. The detrended fluc- 
tuation analysis (♦), the variable bandwidth method (•) and 
the power spectrum density analysis (■). The solid line is the 
roughness exponent for a Levy-flight, C,lf ~ 1/a 

tions we have done roughness exponent measurements 
on random walks on which we impose a power law jump 
distribution with different exponents. The length of the 
samples were 1024 and the number of sample were 1000. 
For power laws with exponents in the range a € {0.5, 3.0} 
we observe corrections for the local window methods and 
the averaged wavelet coefficients method. As seen in Fig. 
Unithe corrections are less than 0.05 for a > 2.0 and in- 
creasing for smaller values of a. A Levy-flight is defined 
as 

I 

h{l) = ^Ahi (15) 
1=1 

where Ahi are independent increments following the 
power law distribution described above. From Bouchaud 
and Gorges [11] one has that 

z«(Z) = (|/i(;)-/i(/)|2) ocZi/«, (16) 

which gives a roughness exponent (lf = cT^ . When 
a > 1.2 C measured with the local window methods and 
the averaged wavelet coefficient method follow C,hF within 
0.05. For oi < 1.2 ^ can not keep up with the increasing 
Clf, which also increases beyond 1. This is as expected 
for the variable bandwidth and the Max-Min methods 
as they are restricted to measure roughness exponents in 
the range C. — [0, 1]. The power spectrum analysis and 
the second order correlation function however measures 
the correlations in the sign change and will only measure 
^ = 1/2 for any a value as the jumps are uncorrelated. 

As seen above two completely different profiles can 
with some methods give the same roughness exponent 
when there is power law distributed noise in the profiles. 
In addition this power law noise can give multi-affine cor- 
rections to the self-affine scaling. 

Unfortunately fracture profiles from experiments and 
simulations often have noise which partly or completely 
disturb the self-affine property of the profile. This noise 
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FIG. 26: Rk/Rk for random walks of length 1024 with a 
power law step size distribution with exponent a) 1.5 and b 
3.0. 



can come from the discretization of the data during 
recording of the profile, the grain or fiber size of the 
material or from overhangs in the fracture front. For 
numerical fracture models recent results show that the 
overhangs in the fracture front have a power law size 
distribution, [l^, [3 In Barabasi and Stanley ^ , rough- 
ening of a surface with an uncorrelated power law noise 
is shown to give multi-affinity below a correlation length 
Ix which depends on the strength of the power law noise 
given by the exponent a, i.e. p{Ah) cx A/i~", and in 
[ig | Mitchell shows that discontinuities in a self-affine 
surface give multi-affinity. The random walks with a 
power law jump distribution a have multi-affine behavior 
on small scales. The region with corrections to scaling 
due to multi-affinity changes with a as seen in Fig. 1261 
Here one sees that the profiles are multi-affine for scales 
smaller than 20 for a ~ 3.0 and not self-affine at all for 
a = 1.5 
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VI. SUMMARY 

For the system size we have studied here the power 
spectrum analysis and the averaged wavelet coefficient 
methods gave the best estimates for the roughness ex- 
ponent over the range of roughness exponents studied 
here. While the averaged wavelet coefficient method is 
reported to give more accurate results for a smaller num- 
ber of samples, [§] this method is prone to systematic 
errors from power law noise. This also applies to the lo- 
cal window methods. We have also seen that using the 
different methods might give roughness exponents that 



differs with a much as 0.1 with additional errors. In ad- 
dition several of the methods have systematic errors that 
varies with the value of the roughness exponent. 

As previously stated by Schmittbuhl et al. 0] one 
should use several different methods for measuring the 
roughness exponent. Researchers measuring the rough- 
ness exponent should also take great care in uncover the 
noise present in the surfaces that are to be studied. This 
to make sure that the measuring methods chosen are ca- 
pable of measuring the roughness exponent properly. 

Acknowledgments: We want to thank Ingve Simonsen 
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